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Abstract 

Dissipationless shock waves in modulational unstable one-dimensio- 
nal medium are investigated on the simplest example of integrable focus- 
ing nonlinear Schr'odinger (NS) equation. Our approach is based on the 
construction of special exact solution of the Whitham-NS system, which 
"partially saturates" the modulational instability. 

1 Introduction 

It is well known that focusing NS equation 

iu t + u xx + 2 | u | 2 u — 0, (1) 
may be rewritten as hydrodynamic type system 

/t + 2(/-«) B = 0, v t + 2vv x -2f x = (2f xx /f-f2/f 2 )j4, (2) 

where 

u = Vf ■ cxp(iip), v = (p x . (3) 
We consider (l)-(3) with the following step-like initial condition 

rjz — ttt f 2,x<0 f 2a,x<0 , „ 
v7(M) = { ^ Q ~{ fl ;,> ,v(x,0) = 0. (4) 

It is natural to call the solution [f(x, t), v(x, t)] or u(x, t) of the problem (l)-(4) 
the shock wave solution. 

The main purpose of this Letter is to describe the qualitative behavior of 
the shock wave at t ^> 1. We note that this problem is much more difficult than 
the analogous problem for stable models [1,2]. 
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The typical and simplest shock problem in modulational stable model was 
first studied by Gurevich and Pitaevskii [1] on the example of Korteveg-de Vries 
(KdV) equation. In [1] modulational theory developed by Whitham [3] was 
used. It is important that the Whitham system for KdV equation as well as for 
other stable models [3] is hyperbolic, which means that Whitham characteristic 
speeds are everywhere real: Im(Si) = 0. Howere for unstable models these 
speeds are generally speaking complex [3]: Im(Si) ^ (elliptic case). This 
complicates very much any analytical investigation of unstable models due to 
exponential growth of linear perturbation of the (locally) constant background. 
In particular there are obvious obstacles in studying of the shock problems 
in unstable systems, which are closely related to the classical problem of the 
modulational instability of monochromatic wave in (1) (see for example [4]). 

In this Letter we demonstrate how to overcome these difficulties on the con- 
crete example (l)-(4). We shall show that famous strategy [1] may be modified 
and successfully applied to (l)-(4). Our approach is based on the conception 
of "partially saturating modulational instability solutions" of elliptic Whitham- 
NS equations proposed in [5] for description of modulational instability in (1). 
The key role in our construction plays a special subclass of Whitham-NS sys- 
tem solutions with vanishing imaginary parts Im(Si) = of some (not all!) 
characteristic speeds (see (12), (17)). 



2 One-phase NS equation solution and its Whi- 
tham-NS modulations. 

It is well known that one-phase solution of (1), (2) has form 

fi(0) = h + (A - h)dn*{^jr^h ■ 0; m}, m = (/i - / 2 )/(/i - h), (5) 
v I (9) = U/2-A/f(9); 9 = x - Ut, (6) 

where /i > f > fi > > fa, A = \J—f 1/2/3 > 0,dn - is the Jacobi elliptic 
function. Elliptic spectral curve corresponding to (5), (6) has branching points 
A;, i = 1,2,3,4; A 2 = A*, A 4 = A* such that (c.f. [6]) 

X 3 = f3-iS = U/4 + ^h/2-i(^h-^h)/2. {l) 

Whitham-NS equations for (1), (2), (5)-(7) can be presented in diagonal form 
(c.f. [6]) 

d\ t /dt + S t {\)d\ t /dx = 0, i = 1,2,3,4, (8) 

5i = C/ + 2Ai 2 /(l- / iA 32 /A 3 i), 5 3 = C/ + 2A 34 /(1 - M14/A13), rq , 
S 2 = St, St = S^, (j, = E{m)/K{m), W 

where Ay = A^ — A^, K, E are the complete elliptic integral of the first and 
second kind respectively, m = A 2 iA43/A 32 Ai4. 
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3 Solution of Whitham-NS system and long-time 
behavior of the NS-hydrodynamic shock wave. 

We construct an approximate solution of the shock problem (l)-(4) for t S> 1 
as follows. In the "external" region on (x,t) plane: (x < 0,x > x + (t)), the 
solution is zero-phase one: 

l,x>t+{t) ' v °( x >°)=°- ( 10 ) 

In the oscillation region < x < x + {t) the "internal" solution is given by 
the one-phase solution (5), (6) with modulated due to (8), (9) parameters. In 
this region we use special solution of the Whitham-NS system (8), (9): 

Ai = const., (11) 

/m(5 3 )=0, t-Re{S 3 )-x = g{(3,5), (12) 

where g((3, 5) is an arbitrary smooth function of its variables, which is deter- 
mined from the initial conditions. 

Let us consider the simplest case g = (c.f. [1], [5,7]). From the initial 
condition (4) we get 

7=1, a = 0, <^=> Ai = -i. (13) 

Proposition. The system (12), (13) with g = is compatible and has unique 
solution with S > 0, (3 > in the "internal" region < x < x + (t). Near the 
boundary x + = x + (t) the solution of the system (12) has the form 

x + = 4\/2£, x = x + -x', < x' < 1, , . 

0Ril/y/2- 7x'/m, S 2 w x'/2V2t. ( ' 

At the boundary x = x + (t) the solution (//, i>j) from (5), (6) is continuously 
glued with 1 from (4). If (x/t) — » +0 the points (A 3 ,A 4 ) closely come to the 
points (Ai, A 2 ). In this limit our solution (//,«/) from (5), (6) degenerates into 
the stationary soliton (breather), corresponding to A 2 = A 4 = i, Xi = A 3 = —i. 

Remark 1. Due to Galilean invariance of (1) the above analysis may be easily 
extended to the case of the shock wave with more general initial conditions 

\/7(M) = { l x x t°o ' ^°) = 2& - (is) 

The corresponding changes in the formulae (11)-(14) are as follows: x + — > Abt + x + , Ai — >6 + Ai i A 3 — > & + 

Remark 2. The mirror symmetry a; — > —a; allows to solve the shock problem 
symmetric to (4) 

V7M) = { 2^>o - w ^°) = °- ( 16 ) 
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We do not go into the details, but just note that corresponding solution of 
Whitham-NS system (8) has the form 



A 3 = const., Jm(5i) = 0, t ■ Re(Si) - x = 0, (17) 

and the picture of the solution u(x,t) at t 3> 1 is symmetric with respect to 
x — > —x transformation to the above picture. 

Remark 3. The problems (4), (15), (16) are of course the simplest examples 
of the shock problems in unstable model (1), (2). As we have shown above 
in these simplest cases the traditional ideology of [1,2] can be applied (after 
proper modification) . For more general initial data it is necessary to invent new 
technique. For example, how to solve the following shock problem 



u(x, 0) 



l,x < ? 
0,x > ' 



We hope to answer this question in the nearest future. 

Remark 4- Note that the above shock wave process (1)-(14) has the same 
typical time t » 1 as ordinary modulational instability (| Im(Sx) |— 1 as £ = 0). 
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